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POR£UDRD 


This voluM of Ch« "Onboard Utilization of Ground Control Points for 
Imagt Corractlon" final raport provldas an axacutlva ausMiy of tha 
study and sloulatlun rstults and Includas our racownandatlons ralatlng 
to futura mlsafon raqulramants. 

Thraa othar volumas !iava baan Incorporatad Into tha final raport. Vol 
II provldas a datallad dascrlptlon of tha study and simulation rasults 
and Volumas III and IV ara appandlcas containing a dascrlptlon of all 
softvara daslgnad and utlllzad undar this contract. 



II. 


INTRODUCTION 


Several current and anticipated tranda threaten to limit the uaefulneaa 
of future NASA remote eenalng miaalona if we continue to employ current 
data handling methoda. Forecaats of Instrumentation trends show an 
increase In sensor resolution from 80 meters for Landsat 1 to 15 meters 
for the Operational harth Resources Setoillte (OERS) planned for launch 
In 1990 (Ref 1). In addition to the Increased reaolutlon, the OERS 
design currently calls for 20 different bands. This translates to an 
Increase in total data rates for earth sensing missions from 1 million 
bits per second (Mbps) to over 100 Mbps; an Increase of two orders of 
magnitude In less than five years. In addition to affecting the data 
rates, Increased sensor resolution affects the following areas: 

o Greater processing requirements for Image correction, formatting, 
and Information extraction. 

o Increased navagatlon accuracy requirements tc provide the necessary 
image distortion coefficients. 

o Increased archiving requirements to accept the larger data volume. 

Increased resolution of the science sensor is a positive trend for 
applications. However, the acquisition of these data Is totally non- 
determlnlstlc In that nothing is known about the quality, content, or 
location of the imagery prior to or even shortly after it is obtained. 

As a result very little of the information is used by the scientific 
community because of such undesirable effects as cloud coverage, 
unwanted scene content, or exposure dates Chat do not coincide with 
those desired. In fact, less than 1% of all previously acquired re- 
motely sensed data have been processed by the user. 

In addition to the data deluge problem. It is currently not feasible to 
exploit Landsat data for such real-time applications as forest fire 
detection and .nonltorlng and flood detection. In fact, with the 
exception of applications in which the observables do not change 
dramatically with time (such as oil exploration), Landsat ’s usefulness is 
limited from an operational standpoint because the data are already stale 
by the time the user gets them. 

User requirements are another important trend associated with remote 
sensing missions because in the long run these requirements play a 
major role in defining the mission. Current data dissemination 
techniques are limited to a several-month turnaround. However, many 
users require this time to be cut to hours , and a set of future users 
are requesting real-time control of the science instrument. The primary 
limiting factor in data turnaround time is the processing required for 
image registration. The impact of non deterministic data acquisition 
is that excessive processing is required to handle these data and 
register unwanted scenes. Excessive processing leads to both rising support 
costs and unacceptable data turnaround.. 
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Th« final trand chat muat ba ovarcoma to raallsa futura mlaaiona la tha 
rapidly aacalatlng coat of groun<^ support. Tha primary drivara of chia 
trand ara manpowar coata and aophxatlratad aqulpmant. 

A datailad axaminatlon of tha prohlam ravaals chat tha following araaa 
ara tha moat critical limitationa of currant ramota aanaing mlaaiona. 

o Ramota aanaing mlaaiona do not amploy any automatad tachnicuoa of 
data avaluacion for alchar almpl > annotation or data nagaticn. A 
simp la datarmlnacion of parcanc«sa cloud covaraga in a acana could 
affactivaly aliminata on tha ordar of SOX of Che data acquirad. 

o Uaara ara subject to tha achadula of tha mlaalon rathar Chan the 
mission being tailored to tha specific acquisition requlramants of 
tha user. This results in tha acquisition of tremendous volumes of 
unwanted imagery alchar due to poor quality caused by atmospheric 
effects. Che geographic location, or the time of acquisition. Also, 
if a user schedules an experiment around the Landsat schedule and 
weather effects are undesirable for chat period of time, it is 
either necessary to wait eight days for the next Landsat over- 
flight or delete remotely sensed data from the experimental re- 
quirements. Unfortunately this happens all too often. 

o All data, including tracking data from remote stations, must be 
sent to a central facility for image correction, annotation, and 
packetization, and then sent to a second facility for archiving and 
eventual dissemination to the users. This results in an unnecessary 
data link and awkward procedures to ensure critical timing require- 
ments are met. This process is neither cost nor time effective. 

o No advantage is being taken of the onboard navigation capability to 
limit the magnitude of image distortions. By providing onboard 
control, the magnitude of the processing Involved can be dramatically 
reduced. 

The purpose of this contract was to Investigate a new approach to re- 
riiote sensing that would meet future mission requirements by overcoming 
these probl£tU£i. 

The primary sources of image distortion are sensor-peculiar errors, 
viewing perspective, coordinate transformation errors, and spacecraft- 
induced errors. With the development of the multilinear array, the 
primary sensor-caused distortions will be the individual placement of 
the detector elements, optical effects, and orientation of the array 
relative to the sensor coordinate frame prior to flight. These errors 
remain fairly constant over time so the resulting distortions are 
deterministic. In referring to Figure Zl-1, viewing perspective is a 
well-known function of local earth radius and vehicle altitude and can 
be computed as shown. 
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It !• aasunMd that th« local earth radlua la known ao thia diatortior 
la alao datarmlnlatlc. Tha primary arror aourca raoalnlng, tharafora, 

Is apacecraf t-lnducad. Tha apacacraft arror aourcaa can ba catagorlsad 
as : 

1) Attitude datanulnatlont 

a) Star tracker accuracy, 

b) Star tracker configuration, 

c) Knowledge of star tracker mlaallgnmant , 

d) Error In star catalog, 

e) Gyro noise, 

f) Uncertainty In gyro bias, nonorthoganallty and misalignment, 

g) Numerical accuracy; 

2) Ephemerls prediction, 

a) Global positioning system (GPS) accuracy, 

b) Dynamic model accuracy, 

c) GPS update Interval, 

d) Numerical accuracy ; 

3) Transformation error between Inertial and earth-fixed coordinates, 

a) Knowledge of UTI, 

b) Knowledge of earth preccesslon, nutation, polar wander, and 
tidal deformation, 

c) Numerical accuracy; 

A) Misalignment between sensor and body coordinates, 

a) Knowledge of linear array orientation, 

b) Accuracy of thermal distortion model, 

c) Vibration modes between two coordinates, 

d) Calibration technique and frequency, 

e) Numerical accuracy. 

For the sake of discussion, assume that all the error is due simply to 
the attitude determination system. To achieve the temporal registration 
requirement of 15 meters, it will be necessary ;o predict attitude to 
within 4 arc seconds as Illustrated in Figure 1-2. The accuracy of the 
state-of-the-art systems using the NASA standard star tracker and gyro Is 
15 arc-seconds (?a) (Ref 2). 
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Even with the development of charge-coupled device (CCD) star trackers, 
the attitude determination capability will be around 6 arc seconds (2o)> 

Note that 2a nuiobers have been used here corresponding to 95% of the 
data. If Ic numbers corresponding to 67Z of the data are used, the 
accuracy goals can be met. However, by adding even one more error term 
such a& misalignment between sensor and body coordinates of 2 arc seconds 
(2-axls alignment accuracy achievable with optical alignment cubes) , the 
total error budget Is exceeded. From the previous discussion, which Ignored 
many error sources, it Is clear that another approach Is required. 

The key to real-time image correction ll<is In being able to accurately 
determine the location of the science sensor's boreslght In earth-fixed 
coordinates. With this capability It : s not only possible to provide 
real-time measurements of the Image distortions, but with the advent of 
the multilinear array (MLA) it may be possible to provide real-time 
image correction using either resampling or special design of the MLA 
focal plane. This capability also provides the heart of a pointing 
system capable of deterministically acquiring Imagery at specific 
'iarth- fixed coordinates. 

Shortly after definition of the feature identification and location 
experiment (FILE), Martin Marietta began the development of a landmark 
tracker or GCP detector centered around experience gained with terminal 
guidance systems. The prlmar/ purpose of the landmark tracker Is to 
provide periodic measurements of the science sensor's boreslght 
position to be used as input to a navigation system. Previous studies 
(Ref 4 thru 7) have shown that the landmark tracker cannot solve for 
both position and attitude without supplemental measurements from 
another source. Another limitation of the landmark tracker operating 
In the visible spectrum is that observations can be obscured by clouds 
and the correlator will lock onto a false target. For these reasons 
the conceptual navigation system consists of the landmark tracker, a 
GPS receiver, and two NASA standard star traclcers to replace attitude 
measurements when the landmark is obscured. 
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Th« purpose of this contract was to slmulatt operation of a navigation 
system centered around image correction, end analyze performance of the 
system under a variety of conditions. Also of interest was a 
sensitivity analysis to determine the optimal sampling interval for 
each sensor, and the total number of landmarks required to satisfy 
mission requirements. 
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III. ONBOARD IMAGE CORRECTION SYSTEM OVERVIEH 


Th« onboard Imag* corroetlon nyttaa haa baan aaparatad Into thraa prl* 
maty subayacama aa llluatratod in Flgura III-l. Tha navigation proc- 
aaaor accapta nMaaurcmanta from tha GPS racalvar, tha atar trackara, 
and tha landmark trackar» and aolvaa for apacacraft atata. Tha ragla- 
tratlon procaaaor la raaponaibla for corralatlng llva imagarv with 
atorad raf.'artnca grotmd control polnta to produca a landmark alghtlng 
vaetor. Thla procaaaor la alao raaponaibla for providing maaauramanta 
of cha dlatortlon coafflclanta. and for providing Imaga raaampllng if 
thla oparatlon la parformad on board. Tha pointing mount controller 
la uaad to laolata tha aanaor from tna attitude limit cycle and to 
provide additional flexibility through aanaor pointing. 
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Figure III-l Onboard Image Covveation System Block Diagram 
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A. REGISTRATION PROCESSOR 

Flgur* III-2 lUustr«c«a th« functional flow of tha raglattatlon proc- 
aasor. A digital rapraoantatlon of all tha ground control points la 
atorad on a high-danalty tapa racordar (or bubbla manory) and raad into 
random accaaa manory aa raqulrad. Bacauaa rf tha raliablllty conaidar- 
ationa aaaociatad with tha racordar, a nuabar . CCPa will ba raad into 
RAM at ona tima. Aa tha apacacraft procaads in orbit, tha landmark 
that ia to ba vlawad naxt la datarmlnad and the digital rafaranca imaga 
la axtractad from RAM. Information atorad in tha GCP flla, along with 

uncartalnty in vahlcla poaltlon and attituda, halpa to dafina tha loca- 

tion and alza of tha aaarch araa. For an oparational ayatam it ia ad- 
vlaabla to uaa a fixad-aiza aaarch araa. Slnca a atar trackar ia baing 
uned for backup, tha maxlmtnn attituda arror will ba in tha vicinity of 
15 arc aeconda l?o) and tha global poaltioning .. 'atam will produca an 
arror around 12 matara. Tharafora a aaarch araa chat ia 32 pixala 
larger than the GCP alza will provide far more araa Chan raqulrad. For 
thla raaaon, and tha fact that a GCP size of 32 producea tha beat ra- 
aulta, a aaarch area alza of 64 was choaan. 

Proceeding with regiatration procsaalng, a teat ia initiated to datar- 
mlne if Che location of the present scan line is coincidant with the 

pofllcion of the search area. If it ia not, a new satcllita state is 

obtair.ad from the navigation processor. This new state is than trans- 
formed into sensor pointing information for use by the pointing mount 
controller. In Che normal mode of operation, the pointing mount controller 
uses gyro data to compen*ace for the attitude limli cycle thus reducing 
cho magnitude of the di^/c /rtione. If the pointing mo int is baing used 
to provide deterministic data acquisition according to aarth-flxad coordinatea, 
the boreslght position provided by the registration processor will also be 
used as a reference to point from. Following tha determination of boreslght 
position, the loop continues until the current scan line is coincidant with 
the search area. 

At this point, image scan lines extracted from the sensor are stored in 
a buffer and become the search area data. The search area is then 
registered with respect to the GCP and the registration vector passed 
to the navigation processor, which uses this information to update the 
estimate of vehicle state. The entire process is then repeated for the 
next landmark. 
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Hifh-Dan«lcy T«p« Racordar 



Figure III-2 Functional Flou of Registration Processor 
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B. 


NAVIGATION PROCESSOR 


figur* IXX-3 lIlutcraCM eh« functional proctaalng porfoioMd by th« 
navigation procasaor. Tha moat critical function parfoctMO la tha 
aaquanclng of maaauramanta and auparvlaory control of tha raglatratlon 
procaaaor. Tha aaquancar flrat dacldaai baaad on a maaouramant pro- 
file » what type of maaauramant ahould ba mada and whan tha maaauramanc 
ahould ba mada. Tha typaa of maaoMramanti tha aaquanca tablaa aatab- 
llsh ara atar trackara 1 and 2 and GPS. After determining tha time of 
tha next maaauremant« tha aaquancar dacldaa, baaad on tha total number 
of scan Unas left to tha landmark and the scan rata of tha sensor, 
whether there la enough time before tha landmark sighting. Xf thaia 
is, tha roaasuramanf Is prccassad normally. If there Is not, control is 
defaulted to tha raglatratlon processor until tha registration vector 
if returned. 



Figure III- 3 Navigation Proaeaaor Funotional Deaaription 
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Afttr the mcasurviMnt is obtAlnad, th« ••tlnuitad itata vector li lnt;«r- 
gratcd forward to the meaeurement tine. Thle proceaa la part of the 
extended square root navigation filter but la Inplenented In a separate 
package to simplify control. The raw naaaurenenta are processed to 
compensate for knowledge of such error terms as bias and misalignment. 
These compensated measurements are used as Input to the Kalman filter 
package that then estimates the spacecraft state. The spacecraft state 
Includes: 

1) Position; 

2) Velocity; 

3) Attitude; 

4) Gyro drift, nonorthoganallty , and scale factor; 

5) Science sensor misalignment. 


The navigation processor Is responsible for sending the navigation 
state to the registration processor at the beginning of each scan 
line. Thils Information Is v’scd to compute the distortion coefficients 
and to control ihc sensur pointing mount. Attitude Is propagated be*- 
tween measurements by a NASA standard gyro package. The pulse train 
output from the gyros is compensated for knowledge of bias, nonorthogo- 
nality, and scale factor and then Integrated to provide vehicle a':titude. 


C. POINTING MOUNT CONTROLLER AND LINEAR ARRAf 

The pointing mount controller cen be Implemented In two different ways, 
depending on the conf igviratlon of the image correction system. If a 
curved focal plana array is used, tha pointing mount controller would 
be responsible for maintaining the boreaight of the instrument along 
nadir. In a system where image correction is provided through resam- 
pling, the pointing mount controller could be used to provide deter- 
ministic data acquisition by pointing to specified eartn-flxed 
coordinates . 

In a curved focal plane cor figuration, che pointing mount controller 
would form an error vector that represents the angular offset between 
the boreaight of the linear array and nadir. This error vector would 
then be used as an input to the pointing mount to correct for the 
error. In this way the viewing angle of the MLA would be isolated from 
the spacecrafL motion caused by either attitude lim.'.t cycles or atti- 
tude corrections. 

In a selective acquisition mode the sensor boreaight position would be 
compared to the desired pointing angle to generate an offset that would 
be used as an input to the positioning mount. Again the viewing angle of 
the sensor can b<s Isolated from the spacecraft motion to simplify the 
image correction process. 

The idea behind a curved focal plane array originates from the under- 
standing that, with the development of the MLA, the primary image dis- 
tortions will be caused by the viewing angle and the curvature of the 
earth (Fig. II-l). This nonuniform sampling is relatively constant and 
is symmetrical about a vector pointing along nadir. It was therefore 
felt that constructing an array with nonuniform sampling for each ele- 
ment that compensated for these effects could eliminate along scan re- 
sampling, which is one of the largest bottlenecks of remote sensing 
(Figure III-4). 

The limitations of this approach are that the central element must al- 
ways point at nadir to within 0.5 the resolution angle of the sensor. 
While this is antlclp'>.ted as being feasible with either Gimbalflex or 
the angular suspension pointing system, a constraint ia placed on the 
system in that it can no longer be used for deterministic data acqui- 
sition through sensov pointing. On the other hand, the cost savings and 
increased data turnaround time make this an attractive possibi''ity that 
warrants further investigation. The key cons Iderat Iona that need to be 
addressed are: 

1) Can the curvature of the focal plane maintain a precision that en- 
sures additional distortions are not introduced? 

2) Is the variation in the local curvature of earth great enough to 
introduce errors greater than 0.*' pixel? 

3) Are variations in the orbital altitude great enough to cause 
distortions? 
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4) Wiat !■ tha coat aavlngs of performing raal-tlma Imaga corractlon 
using this technique varaus complax raaampling proceaaora. 

5) How raal is the daalre to provide deterministic data acquisition 
through sensor pointing? Are the cost savings of deterministic 
data acquisition greater than the cost of producing the curved 
focal plane array? Can the severe distortions Introduced by 
pointing 30»^ off axis be corrected? 

The answers to these questions will form the basis for selecting the 
system architecture to be used on the next generation of remote sensing 
vehicles. However, even If a curved focal plane is not chosen, the 
separation of the science instrument from the motion of the spacecraft 
through the use of a pointing mount will reduce the distortions and 
simplify the process of image correction while providing deterministic 
data acquisition. 


Sensor 

Coordinate Frame 



Ftj'.ife III- 4 Ciu'vad Fojal Plan qrmy 
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IV 


GCP NAVIGATION SYSTEM SIMULATION 


Thtt GCPSIM program has bean configured to provide acientific aimule- 
tiona CO predict tho performance of the GCP detection ayacem. Two 
typea of aimulation modea have been incorporated. The firat r'tde aimu- 
latea the apacecraft, the apacecraft environment and all meaaurementa . 
The aecond mode extracta ground control point meaaurementa directly 
from Landaat Imagery. 

The program waa aet up to provide extensive error analyaia rather than 
simply a convariance analysis. Although a convarlance analysis provides 
a great deal of information, in many cases this information can be in~ 
accurate or misleading. For example, there are cases where the covari- 
ance matrix converges over a period of time while the actual state es- 
timate diverges from the true state. For this reason the simulatl‘>n 
models the true state of the vehicle using the dynamic equations ot mo- 
tion. By propagating the true vehicle state, one can perform a covari- 
ance analysis and also generate the actual state error and measurement 
residuals. 


A. SIMULATED GCP MEASUREMENT 

GCPSIM was designed to provide the ability to analyze the effect of 
various measurement sequences. This was particularly Important when 
studying the effect of GCP spacing, missed GCP sightings, and the ex- 
pected accuracy after traversing a large body of water. The measure- 
ment sequencer designed for GCPSIM allows any mixture of GCP, global 
positioning system, or star tracker measurements and time delays (pe- 
riods during which no measurements were made) of any length. The se- 
quencer determines the type of measurement and the time at which the 
measurement should be made. The true vehicle position state is then 
propagated forward to this time by integrating the nonlinear equations 
of 'tion with some additional process noise to account for modeling 
errors. The attitude state is propagated by looking up the body rates 
in an attitude profile table and integrating these rates. 

The true vehicle state is used along with a measurement model to gener- 
ate an ideal measurement vector. The ideal measurement is then cor- 
rupted with noise, bias, and misalignment terius and compensated for 
knowledge of these values. This allows a careful analysis of the ef- 
fect of misalignment on the state solution. It is important to under- 
stand the effect of bias and misalignment between the landmark tracker 
and body axis because this is the largest unknown factor contributing 
to a pointing error. It is possible to provide frequent calibration 
for these misalignment errors, but it is difficult to model, for any 
length of time, the various processes that cause the misalignment. For 
example, thermal gradients across the vehicle and vibrational .^odes 
within the fle:.;ible structure are complex functions of such things as 
structural design, sun angle, physical properties of the material and 
many other factors. These processes are the most difficult and least 
understood of all engineering problem... 
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Tht compansatcd mcasuramanca art usad as inpuCa into an axtandad 
aquara root Kalraan filtati which aatimataa tha trua vahicia atata. Tha 
axtandad filtar propagataa tha aatimatad navigation atata, tha atata 
tranaition matrix and tha procaaa noiaa array batwaan maaauramanta by 
intagrating tha varioua diffarantial aquationa uaing a fourth-ordar 
Runga Kutta Gill procaaa. Tha aatimatad attituda atata ia propagatad 
by a gyro modal that corrupts tha output with gyro drift, noiaa, non- 
orthogonality, acala factor and miaalignmant. Tha gyro output ia coro- 
panaatad in a similar fashion to tha maaauramant modal. 

Tha aatimatad atata is usad to form an aatimatad maaauramant, which in 
turn is subtractad from tha true maaauramant to obtain a raaidual.. 

The maaauramant v'aaidual and calculatad Kalman gain ara ubad to update 
tha state estimate, which can be compared with tha trua atata to yield 
tha state error. Tha entire process continues until the spacecraft is 
propagatad forward to the run stop time. 

1. Measurement Models 


The types of measurements modeled in GCFSIM Include GPS position and 
velocity, star tracker sightings, and landmark tracker sightings. The 
design philosophy behind the measurement models is that the actual ve- 
hicle state is used with a geometry model to yield an ideal measurement 
vector. This ideal vector is then corrupted with bias, noise, and mis- 
alignment to provide the actual sensor output. The sensor output is 
compensated for some estimate of the error terms and is then used by 
the filter to estimate the vehicle state. The benefit this design 
approach provides ia that it al’riws a detailed analysis of sensitivity 
to misalignments and compensation i>Mllty. A functional description of 
the measurement models is provided in Figure IV-1. 


Dynamics Model and Navigation State Integrator 

The dynamics model and state Integrator are very closely related in 
that the dynamics model calculates the values of the differential equa- 
tion that is used by the integrator. For this reason, the two are dis- 
cussed together in this section. 

a. Dyp' nics Model - The dynamics model calculates the derivative of 
the spacecraft navigational state, which will be integrated to produce 
the navigational state vector. This is done in part by calculating the 
total acceleration of the spacecraft due to solar pressure and gravita- 
tion effects of the sun, moon, and earth, including fourth zonal har- 
monic terms. 

The sequential process for determing the forces acting on the space- 
craft is: 

1) Calculate Julian data; 

2) Calculate solar and lunar positions; 

3 ) Calculate solar pressure acceleration; 

4) Calculate gravitational acceleration due to sun and moon; 

5) Calculate gravitational acceleration due to earth using four zonal 

harmonic terms; 

6) Calculate total acceleration. 
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b. Navigation St* f Infgrator - Th« position teats la advancad In 
tlnw by nuiBarlcal Intagratlon of tha aquae Iona of motion. Tha aacond- 
ordar aquations of motion ara composad of tha axCamal forcaa acting on 
tlia apacacraft. Tha axtarnal forcaa consist of gaopotantlal, lunar and 
aolar gravitation, and radiation prasaura. Tha Intagratlon routlna la 
uaad to propagata both tha actual atata and tha aatlmatad state betwaan 
maasuramant times. A study aho\»d that tha Runga Kutta Gill (RKG) 
fourth-order numarlcal Intagratlon method Ir optimal for this applica- 
tion. It la aelf-startlng, handles variable step sizes, and la suffi- 
ciently accurate. 

3. The Extended Carlson Square Rcot Filter Mechanization 

The Carlson square root filter that was mechanized in GCPSIM is an ex- 
tension of the conventional Kalman filter where the optimal gain, state 
and covariance updates are determined using a matrix that Is tha square 
root of the conv' tttlonal covariance matrix. The advantages of this 
mechanization are higher accuracy for a given machine precision and a 
guaranteed positive covariance. The extended filter Is necessary due 
to the nonlinear nature of the process. The filter propagages the state, 
state transition, and process noise covariance by Integration of differ- 
ential equations rather than linear propagation over the appropriate 
Intervals. 

The sequence process used in the navigation filter Is: 

1) Integrate state to measurement time 

2) Integrate state transition matrix to measurement time 

3) Integrate process noise to measurement time 

4) Propagate covariance matrix to measurement time using state 
transition matrix 

5) Decompose covariance matrix to its square root form 

6) Generate measurement lesldual 

7) Compute Kalman gain 

8) Update state 

B. EXTRACTED GCP MEASUREMEKTS 

The second mode of operation allowed by GCPSIM provides for a mixture 
of modeled and actual measurements. The GPS and star tracker measure- 
ments are modeled as discussed previously. The GCP measurements, how- 
ever, are extracted from actual Imagery obtained from Landsat. 

The second mode of operation functions much the same way as the first 
mode. The sequencer defines the type and time of the star tracker and 
GPS measurements. The GCP data base containing tbe number of scan 
lines to the next GCP Is consulted to see If the GCP will be encoun- 
tered prior to the measurement; If there Is not enough time, the 
measurement defaults to a GCP sighting. The primary difference between 
the two modes of operation is that prior to exiting from the sequencer, 
an Interrupt Is sent via a parallel data bus to a separate task hosted 
on the PDP 11/70 which simulates the sensor system end the operation of 
the registration processor. 
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V. 


GCP ANALYSIS 


For th« purpose of this study, the spscscrsft orbit wss sssumsd to bs 
circulsr with sn Inclination to ths squstorlsl plans at 80 dagraas. 

For aasa of setup and Interpretation of results, the beginning of each 
simulation run was taken to be at a positive ascending mode with the 
Inertial coordinate frame coincident with the earth's coordinate 
frame. This does not detract from the generality of the results since 
specific scenes on the surface or specific star catalogs were not used 
In the simulation. 

A. ERROR BUDGET 

I • Attitude Reference Unit 

The attitude reference unit for this simulation was assumed to be of 
the NASA standard type (DRIRU II). It consists of three nominally or- 
thogonal gyros whose characteristics are fairly well documented. The 
specifications for these gyros were studied and the terms applicable to 
the assumed system configuration and environment were used. The para- 
meters specified and used are outlined In Table V-I. 

2. Star Trackers 


Each of the two star trackers assumed to be on board the vehicle are of 
the NASA standard type. For the purposes of this simulation, the opti- 
cal axes of the two devices were assumed to be orthogonally oriented In 
a plane containing the vehicle position vector and perpendicular to the 
vehicle velocity vector. Documentation, Including the OADS and existing 
specifications, was studied and the error terms felt to be applicable to 
this configuration and environment *rere used. These parameters are out- 
lined in Table V-2. 
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Tab la V~1 Attituda R»fav«na« Unit Et*rov Oa^ffioitnU 


ParaoMCer 

Value of 
Conven** 
tional 
Unitii la 

Value Uaed, lo 

Source 

Gyro Long-Tann Biaa 
Uneartainicy 

0.017 deg/h 

0.08 X 10**8 red/a 

Spec 

Scale Factor Uncertainty 

90 ppm 

90.0 X 10-6 

Spec 

Nonorttiogonality of Mounting 
Package Aligniaant Uncertainty 
Relative to Moi^nting Surface 

2 arc-a 

10 X 10”6 rad 

Aaaumed 

Each Axia 

10 arc*“a 

50 X 10-6 

Aaaumed 

Gyro Noise Equivalent Angle 

1.3 arc-a 

10 X 10-6 rad/a 

Spec 

Absolute Bias 

2 deg/h 

10 X 10-6 rad/a 

Spec 

Absolute Scald Factor 

0.5X 

0.005 

Spec 


Table V-2 Stzx* Tixiakev Ervoi* Coeffioienta 


Parameter 

Value of 
Conventional 
Unite, 1 a 

Value Uaed, la 

Source 

Biaa, Each Axia 

10 arc-*a 

50 X 10-6 rad 

Spec 

Quantisation 

1 arc -a 

5 X 10-6 rad 

0A08 

Misalignment 
Bach Axia 

5 arc -a 

25 X 10-6 rad 

OADS 

Scar Catalog 

1 arc -a 

5 X 10-6 rad 

OADS 


3. Global PjaitioninK System (GPS) 

The GPS is assumed to be available Co the spacecratt virtually continu 
ously. The update interval using GPS was predominantly 3 seconds; how 
ever, variations have been made in this parameter to study its affect 
on system performance. The GPS system error parameter assumed for the 
purposes of this simulation are presented in Table V-3 . 


Table /-J Evvov Caeffiaianta 


Parameter 

Value of 
Conventional 
Unite, In 

Value Used, lo 

Source 

Velocity Biaa 

0 

0 

Assumed 

Velocity Noise 

0.006 m/s 

6 X 10 -6 tan/s 

Assumed 

Position Bias 

0.04 m 

4 X 10-5 km 

Assumed 

Position Noise 

4.0 m 

4 X 10“3 km 

Assumed 
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Landmark TrackT (LMT ) 


Tha LMT waa asaumad to hava It a optical axla collnaar with tha nadir 
pointing axia tha apacacraft. Tha look angle to tha landmark ia in 
the plane containing the optical axia, perpendicular to tha vehicle 
velocity vector. It haa a maximum value of 10 dagreea with two modaa 
of operation. Theac modaa are uaer-aelectabla to take on either fixed 
specified look anglea or uniformly diatributed look anglea of a aped- 
fled maximum magnitude. Table V-A outlinea the error characteriatlca 
assumed foi the LMT. 

Table V~4 Landmark Traakar Exror Coeffiai9nta 


Parameter 

Value 

Conventional 
Units 1 a 

Value Used 1 o 

Source 

Biaa, Each Axia 

2 arc -a 

10 x 10*~6 rad 

Assumed 

Quantization 
Each Axia 

1 arc -a 

5 X 10“6 rad 

Assumed 

Hiaalignment 
(three axaa) 

5 arc-a 

25 X 10"6 rad 

Aaaumad 

Landmark Location 
Uncertainty 

1 arc**8 

5 X 10“6 rad 

Assumed 


PARAMETRIC STUDIES 

The intent of the parametric studies was to determine the system's 
sensitivity to such parameters as LMT update interval, look angle mode 
and range, GPS update interval and accuracy, and recovery time after a 
long period without LMT Updates with and without star tracker backup. 
The study waa progressive in the above order using as parameters the 
values that appeared to be optimum from the previous study. With the 
exception of the study of system sensitivity to the GPS accuracy, the 
error budget referred to was used uniformly for all studies. 

Table V-5 outlines the cases studied by title and most significant 
parameters . 




Tahl« I'-f Studied 



GCP Update 

CPS Update 


TicU 

Interval, a 

Interval, a 

Raiaarka 


GCP Updae« Interval Saneitivlty 


CASE 1 

5 

5 

CASE 1 through CASE 7 ware run 

CASE 2 

10 

5 

with ♦lO-deg random look anglea 

CASE 3 

20 

5 

aa well aa a 0-dogree fixed look 

CASE U 

40 

5 

angle. 

CASE 3 

80 

5 


CASE 6 

160 

5 


CASE 7 

320 

5 


Look Angle Senaitivity 

CASE 9 

80 

5 

0-deg fixed look angle. 

CASE 10 

80 

3 

2.5-deg fixed look angle. 

CASE 12 

80 

3 

3-deg fixed look angle. 

CASE 14 

80 

3 

10**deg fixed look angle. 

CASE 17 

80 

3 

-»'2> 3-deg random look angle. 

CASE IS 

80 

5 

^3-deg random look angle. 

GPS Senaicivit) Update Interval 

CASE 21 


10 

CASE 21 through CASE 26 were 

CASE 22 


20 

run with ♦lO-deg 1 ok angles. 

CASE 23 


40 


CASE 2^ 

80 

80 


Accuracy 

CASE 23 

80 

3 

GPS velocity error increased by 4X. 

CASE 26 

80 

3 

GPS position error increased by 4X. 

Recovery 

Sensit ivity 



CASE 30 

80 

3 

♦10-deg random look angles; linear 
cloud table; default to star 




t rackets . 

CASE 31 

80 

3 

♦lO-deg random look anglea; GCP 
shut down for lOUO s with default 
to star crackers. 

CASE 32 

80 

5 

♦10-deg random look angle; GCP 
shut down for 1000 s without 
default to star crackers. 























c. 


SUMMARY OF RESULTS 


Tha ayatam aa aimulaead OMt the raqulramanra aatabllahad aarllar for con> 
dltlona whara tho GCP updata Intarvala ara batwaan 20 and 80 aaconda 
with GPS updataa takan at about 10-aacond Intarvala. In aoma caaaa part 
of tho initial tranaiant waa includad in aatlmation of tha paak or value 
for a givan paraoator. Tharofora laoat of tha roaulta ara fait to ba 
paaaimiatic with roopact to tho aiza of tha arror boxmda. Thooa roaulta 
do not roflact tha aatlmation of auch gyro parwmotara aa gyro biaa, 
miaalignmant , acala factor arror and nonorthogonality. With oatimatlon 
and companaation for thaoa paramatara, tha roaulta ahould ba furthor im- 
provod . 

Thaae atudlaa took advantage of tha fact that tho rofaranco body 
coordinato ayatam ia definad within tho landmark trackor or GCP 
aanaor. Thorafora tha miaalisnmant of that aanaor ia by definition 
equal to zaro. Howavar, thia ahlfta tha arror contribution of mia> 
alignment of other aanaora if auch mlaalignmanta arc not aatlmatad and 
compenaated. Tha off acta of auch uncompanaatad torma aa gyro biaa may 
bs aeon in caaea where long GCP update intarvala ara uaad. Further 
atudieo ahould be undertaken to implement the aatlmation and compenaa'- 
tion of conatant gyro parametera and poaalbly tha landmark tracker 
parameters of bias and misalignment ahould that sensor coordinate frame 
not be taken as a body reference. 
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VI. CONCLUSION AND RECOWEWDATION8 


A. CONCLUSIONS 

On« of th« koy bottlonocks •avoclotod with NASA's tiu!>to-«nd data 
managamant problam la tha procaaa of Imago corractlon. Raaanpllng not 
only raqulrao a tromandoua amount of data procaaalng but alao roqulroa 
tha intagratlon of volumoa of data from ovary aapact of tha mission, 
Including gyro data, tracking data, and Imagary. 

With tha davolopmant of tha MLA, autonomoua raal-tlma Imago corractlon 
will bacoma faaslblo. Tha distortions assoclatad with tha 
multlspactral acannar's nonllnaar mirror swaap and nonuniform swoop 
parlod will ba allmlnatad, laavlng a wall-daflnad dlotortlon causad by 
aarth curvatura and look angla. Tha tromandoua Incraasas In 
computational capabllltias both on tha ground and aboard tha spacacraft 
will alao simplify tha problam. 

A solution to tha problam may ba roallzod through olthar ground 
procaaalng or spacoborna procaaalng. Howovar, avan though ground-basad 
procaaalng may allovlata tha Immadlato problam. In tha long rango 
adaptlva syatoms using onboard Intalllganco will bo roqulrad. 

Tharefore It may be advantagaous to begin tha transition to onboard 
automated systams now. 

In the past there has been a lack of coordination between the 
scientific user conmunlty and the engineers responsible for spacecraft 
design. This has resulted In a physical separation between tha design 
and Implementation of the science payload and the control system. An 
example of this thinking Is shown In the multlmlsslon spacacraft (MMS) 
where subsystems are treated as modules and the payload Itself is 
physically separated from tha control portion of the vehicle. This 
type of design, although attractive from a standardization point of 
view, Ignores the inherent relationships between uaer requirements, 
spacecraft control requirements, and ground support requirements. It 
is possible that rather than saving significant costa, standardization 
may result In higher end-to-end costs. This design philosophy must be 
reevaluated with regard to future missions. In the end-to-end design 
of remote sensing spacecraft, the primary emphasis In the guidance and 
control system must shift from simply estimating the ephemerls and 
attitude of the spacecraft to estimating the position of the science 
sensor's FOV on the earth’s surface. This provides the basis for both 
real-time Image correction and deterministic data acquisition through 
sensor pointing. This shift of emphasis wll] Impact the design of the 
entire spacecraft. For example. If the science sensor boreslght 
position is to be determined. It Is desirable to place the gyro package 
close to that sensor to reduce the misalignment between the two. This 
(MMS) configuration, which provides a physical separation between the 
payload and the guidance and control system, may not satisfy the 
requirements of many future remote sensing missions. 

Many benefits may be realized through onboard spacecraft control. By 
isolating the science sensor from the attitude limit cycle of the 
vehicle using a pointing mount, a curved focal plane MLA that provides 
uniform sampling of the image can be incorporated, thus reducing the 
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•tnount of rosaspling roqulrtd. This «pproach, alchough attractlva for 
mlaalona auch aa Landaae vhara eha raviait eyela la parlodlc, may not 
ba attracclva for mlaalona raqulrlng datarmlnlatlc data a'^qulalelon 
through aanaor pointing bacauat a curvad focal plana array boraalght 
muat alwaya point at nadir. Howavar, a pointing mount la uaaful In 
althar caaa to laolata tha aclanca aanaor from tha atcltuda limit cycla. 

Ona of tha kaya to futura ramota aanalng mloalona la tha onboard 
datanalnatlon of aanaor boraalght poaltlon In aarth-flxad coordlnataa. 

Be t Imaga corractlon and datarmlnlatlc data acqulalclon ayatama will 
raqulra thla knovladga to ba uaaful. Analyala haa ahown that Inartlal 
ayatoma alona do not provlda aufflclant accuracy for aavaral raaaona: 

1) It la difficult to raaolva mlaallgnmanta batvaan Inartlal aanaor 
framaa and aclan** i aanaor framaa; 

2) Mlaallgnmant batvaan Inartlal and local aarth framaa ara difficult 
to coffiputa: 

3) Star trackara ara not accurata anough. 

Ground control point (GCP) alghtlnga can ba uaad to directly aolva for 
thaoa mlaallgnmanta and hanca reduce tha error In boraalght 
determination. 

Analysis has shown that a ragiatratlon accuracy of 0.1 plxela can be 
achieved ualng an SSDA correlation algorithm. It la alao poaalbla to 
detect altuatlona where the correlator locka onto a falae target. 

Theae factora make the correlator an attractive navigation aanaor for 
an onboard Image correction system. 

A system comprising two NASA standard star trackers, a NASA standard 
gyro package, a GPS receiver, and the landmark tracker with a SO-meter 
resolution was simulated. Analysis results indicate that with a GPS 
update frequency of between 5 and 10 seconds and a GCP sighting 
frequency of between 20 and 80 seconds, the position of the sensor 
bo.’eslght can be determined to within 15 meters. With the one 
exception of the MLA, tha entire system can be implemented with 
existing technology. Tha most restricting feature of the system Is the 
tremendous storage requirements for the GCP data base. Assuming that 
GCP *8 may only be located land (accounting for 20% of the earth's 
surface) and assuming a GC update frequency every 40 seconds or 
approximately every 11/3 scenes, the total storage requirements may be 
found by; 

3 2 2 

Total earth surface area ■ 4-(6 x 10 ) km 

g 

* A • 5 X LO * 3 2 2 2 

Total number of available scenes ■ 4.5 x 10® km^ x 0.2 175 km^ 

- 2.9 X 10^, 

Total number of GCPs ■ 2.9 x 10’’ x 3/4 

- 2 x.^10^’ 2 6 

Total memory required ■ 2 x 10' GCPs x 32 pixels x 8 bits ■ 18 x 10 bits. 

GCP pixel 
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Although this Is s trensndous storsgs rsqulrsmsnc , It is still within 
the limits of current technology This volume msy also be reduced by a 
factor of two by reducing the number of gray scales used in the reference 
GCP from 256 to 16. If histogram equelizatisu is used to reduce the 
number w5 gray scales, the accuracy of the correlation may not be 
affected. This possibility should be investigated. 

It is estimated that processing requirements can be met with the use of 
two general-purpose flight computers and a hardwired correlator. The 
separation cf functions would be directly analogous to the navigation 
and the registration processors discussed in Section II. 


B. RECOMMENDATIONS 


If the landmark tracker is to be realized as an operational system, a 
well-planned development schedule is required. The critical tasks that 
must be performed include: 

1) Further analysis of the correlator is required to determine its 
sensitivity to seasonal variations. It is also necessary to determine 
the types of landmarks that produce the best correlation. This in- 
formation will be invaluable in the selection of GCPs for an opera- 
tional system; 

2) Research must be conducted to determine methods of reducing the 
storage requirements of the GCP data base. The most promising 
approach seems to be a mixture of image enhancement and gray scale 
reduction; 

3) A study must be conducted to establish the configuration of the 
onboard processing network; 

4) The entire system should be breadboarded and ground-tested to 
resolve peculiarities of the design; 

5) The breadboard system may be flown on a low-cost aircraft flight 
test to more fully demonstrate feasibility; 

6) The breadboard system could also be Incorporated into the IAS 
ground demonstration as an additional step toward a total system 
configuration. 

The curved focal plane MLA is another concept that should be 
investigated. If feasible, the concept would dramatically reduce the 
processing requirements associated with remote sensing missions. This 
Investigation can be separated into two tasks. First, the feasibility 
of the concept must be analyzed. Technical considerations include; 

1) Precision required in machining of the curved surface; 

2) Effect of variations in local earth curvature on registration 
accuracy ; 

3) Pointing requirements needed to ensure 1/2 pixel registration. 
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Th« •«cond taak that muat ba parforaad la a coat/bar.af It tradaoff 
batwaan datarmlnlatlc data acqulsltlQn through tanaor pointing varaua 
Imaga corraction using tha curvad focal plana array. Tha two 
approachaa ara mutually axclualva. 
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